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The origin of cosmic magnetism is an issue of fundamental importance in astrophysics. We review 
here some of the ideas of how large scale magnetic fields in the universe, particularly in galaxies 
and galaxy clusters could arise. The popular paradigm involves the generation of a seed magnetic 
field followed by turbulent dynamo amplification of the seed field. We first outline various seed 
field generation mechanisms including Biermann batteries. These in general give a field much 
smaller than the observed field and so they require further amplification by dynamo action. The 
basic idea behind fluctuation dynamos, as applied to cluster magnetism and the mean-field helical 
dynamo as applied to disk galaxies, are outlined. Major difficulties with the dynamo paradigm 
are considered. It is particularly important to understand the nonlinear saturation of dynamos, and 
whether the fields produced are coherent enough on large-scales to explain the observed fields in 
galaxies and clusters. At the same time the alternative possibility of a primordial field lacks firm 
theoretical support but can have very interesting observational consequences. 
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1. Introduction 

Magnetic fields are crucial for understanding a number of physical processes in the universe. 
They give "fluid" like properties to hot and dilute plasma and affect transport properties like ther- 
mal conduction, viscosity and resistivity of plasmas. Present day star formation is controlled by 
magnetic fields with the field possibly determining the typical masses of stars and the initial mass 
function. They are important for accretion and ejection flows from protostars to active galactic 
nuclei (AGN). The galactic field is important for cosmic ray confinement and transport. Their most 
important role for radio astronomers is that they make it possible to "see" objects due to synchrotron 
and inverse Compton emission of relativistic electrons gyrating in the magnetic field. Indeed it is 
this radio emission that is crucial to measure the strengths and coherence scales of magnetic fields 
on the largest scales in the universe, that associated with galaxies and galaxy clusters [JTJ] _ Such 
studies show average total field strengths in nearby spiral galaxies B ~ 10//G with a mean field 
about half this value, correlated on scales of several kpc to tens of kpc. Galaxy clusters also seem 
to host magnetic fields of several /jG, and correlated on several kpc scales. How do such ordered 
large-scale fields arise in galaxies and clusters? Indeed the origin cosmic magnetism in general 
is one of the fascinating questions of modern astrophysics. It forms one of the key goals of the 
Square Kilometre Array (SKA). We review here some of the ideas that have been put forth on how 
the universe got magnetized. 

We start with a brief outline of the basics of magnetohydrodynamics, consider then both seed 
field generation mechanisms and dynamo amplification processes. We consider the possibility of 
magnetizing the universe through outflows from starbursting galaxies and AGN. We finally discuss 
the possible role of primordial magnetic fields created perhaps in an early universe phase transition. 
We focus on magnetic fields on scales larger than galaxies here, as such fields will also be the focus 
of SKA type telescopes, and do not discuss stellar or planetary magnetism. An extensive review of 
several aspects of astrophysical magnetic fields and nonlinear dynamos can be found in 



2. Basic MHD 



The evolution of the magnetic field B is governed by the induction equation 

— = Vx(UxB-nVxB) (2.1) 
at 

got from combining Maxwell equations (neglecting the displacement current) and a simple form of 
Ohms law. Here U is the fluid velocity and 77 the microscopic resistivity. The field back reacts on 
U via the Lorentz force J x B 



P 



dU 

— + (U-V)U 

ot 



= -V/> + — + f + F vlsc , (2.2) 



where p, p are the fluid density and pressure, f represents all the body forces and F v j sc the viscous 
force. A few simple consequences of the above equations are as follows: (i) If U = 0, then the 
field simply decays due to finite resistivity, (ii) On the other hand for rj — > 0, the magnetic flux 
<D = B • dS through any surface moving with the fluid is 'frozen' in the sense that d$>jdt — > 0. 
(hi) The magnetic Reynolds number R m = (UB)/(tjB/L) = UL/rj measures the relative importance 
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of the induction and resistive effects. (Here U and L are respectively the typical velocity and 
length scale in the system.) For most astrophysical systems R m » 1. In galaxies we expect naively 
R m ~ 10 18 and in clusters R m ~ 10 29 based on the Spitzer resistivity Plasma instabilities in the 
presence of a weak magnetic field could produce fluctuating magnetic fields which scatter charge 
particles and affect various transport processes [g]. But to affect resistivity by scattering electrons, 
these fluctuations need to have significant power down to the electron gyro radius, (iv) Note that 
B = is perfectly valid solution to the induction equation! So magnetic fields can only arise from 
a zero field initial condition, if the usual form of Ohm's law is violated by a "battery term" (see 
below), to generate a seed magnetic field B seec i- (v) The resulting seed magnetic field is generally 
very much smaller than observed fields in galaxies and clusters and so one needs the motions (U) 
to act as a dynamo and amplify the field further. It turns out to be then crucial to understand how 
dynamos work and saturate. 

At this stage it is important to clarify the following: it is often mistakenly assumed that if one 
has a pre-existing magnetic field in a highly conducting medium, and the resistive decay time is 
longer than say a Hubble time, then one does not need any mechanism for maintaining such a field. 
This is not true in general; because given a tangled field the Lorentz forces would drive motions 
in the fluid. These would either dissipate due to viscous forces if such forces were important (at 
small values of the fluid Reynolds number Re = UL/v, where v is the fluid viscosity) or if viscosity 
is small, drive decaying MHD turbulence, with a cascade of energy to smaller and smaller scales 
and eventual dissipation on the dynamical timescales associated with the motions. This timescale 
can be much smaller than the age of the system. For example if clusters host a few /uG magnetic 
fields say tangled on say / ~ 10 kpc scales, such fields could typically decay on the Alfven crossing 
timescale I/Va ~ 10 8 yr, where we have taken a typical Alfven velocity Va ~ 100 km s . Although 
the energy density in MHD turbulence decays with time as a power law, this time scale is still much 
shorter than the typical age of a cluster, which is thought to be several billion years. Similarly, if 
the fluid is already turbulent, the associated larger turbulent resistivity will can lead to the decay 
of large-scale fields (except perhaps in the presence of strong shear, where there may be dynamo 
action; see below). Therefore, one has to provide explicit explanation of the origin and persistence 
of cosmic magnetic fields; reference to the low Ohmic resistivity of the plasma is not sufficient 
if the gas is turbulent or the magnetic field is tangled. Hence even if the medium were highly 
conducting, in most cases, a dynamo is needed to maintain the observed magnetic field. 

3. The first seed fields 

We concentrate on astrophysical batteries here and discuss primordial magnetic fields arising 
from the early universe in Section 7. The basic problem any battery has to address is how to produce 
finite currents from zero currents? Most astrophysical mechanisms use the fact that positively and 
negatively charged particles in a charge-neutral universe, do not have identical properties. For 
example, if one considered a gas of ionized hydrogen, then the electrons have a much smaller mass 
compared to protons. Thus that for a given pressure gradient of the gas the electrons tend to be 
accelerated much more than the ions. This leads in general to an electric field, which couples back 
positive and negative charges. If such a thermally generated electric field has a curl, then from 
Faraday's law of induction a magnetic field can grow. The resulting battery effect, known as the 
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Biermann battery, was first proposed as a mechanism for the thermal generation of stellar magnetic 
fields [@]. 

The thermally generated electric field is given by E^ !f , r = -Vp e /en e got by balancing the forces 
on the electrons, due to pressure gradient and the electric field and assuming the protons are much 
more massive than the electrons. The curl of this term leads to an extra source term in the induction 
equation, which if we adopt p e = n e k^I ', where k# is the Boltzmann constant, gets modified to, 

? = Vx(UxB-,VxB)-^xVr. (3.1) 

ot e n e 

Therefore over and above the usual flux freezing and diffusion terms we have a source term which 
is nonzero if and only if the density and temperature gradients, Vn e and V7\ are not parallel to each 
other. In the cosmological context, such non-parallel density and temperature gradients can arise 
in a number of ways. For example, in cosmic ionization fronts produced when the first ultraviolet 
photon sources, like starbursting galaxies and quasars, turn on to ionize the intergalactic medium 
(IGM), the temperature gradient is normal to the front. However, a component to the density 
gradient can arise in a different direction, if the ionization front is sweeping across arbitrarily laid 
down density fluctuations, which will later collapse to form galaxies and clusters. Such density 
fluctuations, will in general have no correlation to the source of the ionizing photons, resulting in 
a thermally generated electric field which has a curl, and magnetic fields correlated on galactic 
scales can grow. After compression during galaxy formation, they turn out to have a strength 
fi~3xl(T 20 G This scenario has in fact been confirmed in detailed numerical simulations 
of IGM reionization [g]. The Biermann battery has also been shown to generate both vorticity 
and magnetic fields in oblique cosmological shocks which arise during cosmological structure 
formation [|7|]. 

The asymmetry in the mass of the positive and negative charges can also lead to battery effects 
during the interaction of radiation with ionized plasma. Note that the Thomson cross section for 
the scattering of photons with charged particles depend inversely on the mass of the particle. So the 
electron component of an ionized plasma is more strongly coupled with radiation than the proton 
component. Suppose one has a rotating fluid element in the presence of a radiation bath. The 
interaction with photons will brake the velocity of the electron component faster than the proton 
component and set up a relative drift and hence lead to magnetic field generation [§]. In the modern 
context, second order effects during recombination also leads to both vorticity and magnetic field 
generation due to the y-e/p scattering asymmetry. The resulting magnetic fields are again very 
small B ~ 10~ 30 G on Mpc scales upto B ~ 1(T 21 G at parsec scales 

As can be gleaned from the few examples considered above all battery mechanisms give only 
very small fields on the cosmological scales much smaller than the observed fields. One therefore 
needs some form of dynamo action to amplify these fields further. Of course, seed fields for 
dynamos in larger scale objects like galaxies and clusters can arise not necessarily due to battery 
effects, but also due to more rapid dynamo generation in objects with a much shorter dynamical 
timescale, like stars and AGN, and subsequent ejection of these fields in to the interstellar and intra 
cluster media (cf. [[ll]] for reviews). In this case much larger seed fields say for the galactic dynamo 
B ~ 10~ 9 G are possible [|TTj] . The down side is the as yet unresolved question as to how magnetized 
gas say ejected in a supernovae or AGN is mixed with unmagnetized gas in the protogalaxy, and 
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how this mixing affects the coherence scale of the field? And even in this case one still needs 
dynamos to work efficiently in stars and AGN to generate the seed field, and in galaxies and clusters 
to further amplify and maintain it against decay. 



4. Fluctuation dynamos in galaxies and clusters 

Magnetic fields in a conducting medium can be amplified by the inductive effects associated 
with the motions of the medium. In this process, generally referred to as a dynamo, the kinetic 
energy associated with the motions is tapped to amplify magnetic energy. The plasmas in galaxies 
and clusters are most often turbulent. The dynamo in this context is referred to as a 'turbulent dy- 
namo', and its analysis must rely on statistical methods or direct numerical simulations. Turbulent 
dynamos are conveniently divided into fluctuation (or small-scale) and mean-field (or large-scale) 
dynamos. The fluctuation dynamo produces magnetic fields that are correlated only on scales of 
the order of or smaller than the energy-carrying scale of the random motions. We will discuss this 
here and turn to mean-field dynamos in the next section. 

The importance of fluctuation dynamos in cosmic objects obtains because they are generic in 
any random flow where R m exceeds a modest critical value R mfiT ~ 100. Fluid particles in such a 
flow randomly walk away from each other. A magnetic field line frozen into such a fluid will then 
be extended by the random stretching (if R m is large enough). Consider a small segment of a thin 
flux tube of length / and cross-section A, and magnetic field strength B, in a highly conducting fluid. 
Then, as the fluid moves about, conservation of flux implies BA is constant, and conservation of 
mass implies pAl is constant, where p is the local density. So B oc pi. For a nearly incompressible 
fluid, or a flow with small changes in p, one will obtain B oc /. Any random shearing motion 
which increases / will also amplify B; an increase in / leading to a decrease in A (because of 
incompressibility) and hence an increase in B (due to flux freezing). 

Of course since the scale of individual field structures decreases, (that is since A ~ l/B), as 
the field strength increases, the rate of Ohmic dissipation increases until it compensates the effect 
of random stretching. For a single scale random flow, this happens when V0//0 ~ 77//^, where vo 

is the typical velocity variation on scale Iq and is the scale of the magnetic field. This gives 

1 11 

Ib- Irt ~ fo/Rm > the resistive scale l v for the flow, where R m - vqIo/J]- What happens after this 
can only be addressed by a quantitative calculation. For a random flow which is delta-correlated in 



time, it was shown by Kazantsev [12], that magnetic field can grow provided R m > /? m , cr ~ 30- 100, 
depending on the form of the velocity correlation function. In the kinematic regime the field grows 
exponentially roughly on the eddy turnover time Io/vq and is also predicted to be intermittent, i.e., 
concentrated into structures whose size, in at least one dimension, is as small as the resistive scale 
Ij] (e.g., [11, |IJ]). In Kolmogorov turbulence, where the velocity variations on a scale / is v/ oc /^ 3 , 



the ^-folding time is shorter at smaller scales, l/vi oc Z 2 ^ 3 , and so smaller eddies amplify the field 
faster. If P m = R m /Re » 1, as relevant for galactic and cluster plasma, even the viscous scale eddies 
can exponentially amplify the field. 

The fluctuation dynamo has since been convincingly confirmed beyond the limits of Kazant- 
sev 's theory by numerical simulations of forced and flows [[14, 15], especially when R m > Re. Such 



simulations are also able to follow the fluctuation dynamo into the non-linear regime where the 
Lorentz forces becomes strong enough to affect the flow as to saturate the growth of magnetic field. 
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In the context of galaxy clusters turbulence would mainly be driven by the continuous and 
ongoing merging of subclusters to form larger and larger mass objects. One typically expects the 
largest turbulent scales of Iq ~ lOOkpc and turbulent velocity vo ~ 300 km s , leading to a growth 
time To ~ lo/vo ~ 3 x 10 8 yr; thus for a cluster lifetime of a few Gyr, one could then have signif- 
icant amplification by the fluctuation dynamo (cf. [|l7|]). And in the case of galactic interstellar 
turbulence driven by supernovae, if we adopt values corresponding to the local ISM of Iq ~ lOOpc, 
vo ~ 10km s" 1 one gets to ~ 10 7 yr. Again the fluctuation dynamo would rapidly grow the magnetic 
field even for very young high redshift protogalaxies. 

The major uncertainty, in case we want to use the fluctuation dynamo to explain observed 
Faraday rotation measures in a galaxy clusters or that inferred in some high redshift protogalaxies, 
is how intermittent the field remains when the dynamo saturates. A simple model exploring an 
ambipolar drift type nonlinearity [16], suggests that the smallest scale of the magnetic structures 

— 1/2 

will be renormalized in the saturated state to become l R - k)R m ,cr instead of the resistive scale l v . 
This essentially happens via a renormalization of the effective magnetic diffusivity in these models 
[16], such that the dynamo is saturated via a reduction of the effective magnetic Reynolds number 
down to its critical value for the dynamo action. In this case one could indeed obtain significant 
Faraday rotation measure through such a fluctuation dynamo generated field. On the other hand 
it has been argued that the fluctuation dynamo generated field remains highly intermittent even at 
saturation, with field reversals typically occurring on the resistive scale [15]. Since the cluster 
R m ~ 10 29 if one uses naively the collisional Spitzer resistivity, one hardly expects to see any 
Faraday rotation from such a field. Plasma effects would then important to renormalize the effective 
R m , much below this ridiculously large value [Q]. One needs a better theoretical understanding of 
the non-linear saturation of fluctuation dynamos and the extent to which plasma effects alter the 
resistivity of cluster plasma, to make further progress. 

Meanwhile direct simulations of the fluctuation dynamo at modest values of R m ~ 400- 1000, 
to examine the resulting Faraday rotation, have been carried out [17], following the methods of 
Haugen, Brandenburg and Dobler [jjjp. The simulations use a Cartesian box on a cubic grid with 
256 3 mesh points. The computational box contained just a few turbulent cells on the forcing scales. 
The Faraday rotation measure K Jn e B-dl for 256 2 lines of sight through the computational domain 
was calculated, and various properties of its distribution were determined including the normalized 
standard deviation cr RM , normalized by Kii e B Tm Jo (Here n e and h e are the actual and average elec- 
tron densities respectively, B rms the RMS value of the magnetic field in the box, lo the forcing scale 
of the random flow, and K is a constant of proportionality which converts the line of sight integral 
into physical units of RM). For a random field one expects a distribution of Faraday rotation with 
mean zero and a standard deviation cr RM , which depends on the degree of the field intermittency. If 

— 1/2 

the field saturated with many resistive scale reversals one expects ctrm 1, whereas if - loR m ,cr 
decides the coherence scale of the field, one expects cr RM of order unity. The simulations found the 
normalized cr RM ~ 0.5 for P m = l,Re = R m « 420 and cr RM ~ 0.3 for a P m = 30, Re « 44, R m ~ 1300 
case. Therefore the fluctuation dynamo does indeed seem to be able to produce magnetic fields 
which will give significant Faraday rotation measure, although simulations at much higher values 
of /? m ,Re are required. Further, as emphasized by [[17]], during epochs when cluster turbulence 
decays, one can still get a significant ct R m, as it decreases in time as a shallower power law than the 
field itself, because of an increase in the field coherence scale. Similar ideas involving cluster tur- 
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bulence and fluctuation dynamos have also been invoked by [18] to explain magnetic fields in cool 
cores of galaxy clusters. The way the fluctuation dynamo saturates will be important not only to 
decide if it can result in observable Faraday rotation, but also to issues of Cosmic ray confinement 
in the early protogalaxy. 

5. Mean field dynamos and galactic magnetism 

A remarkable change in the turbulent dynamo action occurs if the turbulence is helical. This 



can be clearly seen for example in the simulations by Brandenburg [19], where a large scale field, 
on the scale of the box develops when a helical forcing is employed, even though the forcing itself 
is on a scale about l/5th the size of the box. The large scale field however in these closed box 
simulations develops only on the long resistive timescales. It is important to understand how such 
a field develops and how one can generate a large scale field on a faster timescale. The possible 



importance of helical turbulence for large-scale field generation was proposed by Parker [20], and 



is in fact discussed in text books [21]. We summarize briefly below the theory of the mean-field 
dynamo (MFD) as applied to magnetic field generation in disk galaxies, turn to several potential 
problems that have been recently highlighted and their possible resolution. 

Suppose the velocity field is split into the sum of a mean, large-scale velocity U and a turbulent, 
stochastic velocity u. The induction equation becomes a stochastic partial differential equation. Let 
us split the magnetic field also as B = B + b, into a mean field B = (B) and a fluctuating component 
b. Here the average <>, is defined either as a spatial average over scales larger than the turbulent 
eddy scales (but smaller than the system size) or as an ensemble average. Taking the average of the 



induction equation [2.1| , one gets the mean-field dynamo equation for B, 



3D 

— = Vx(UxB + £-??VxB). (5.1) 
at v 1 



This averaged equation now has a new term, the mean electromotive force (emf) £ = uxb, which 
crucially depends on the statistical properties of the small-scale velocity and magnetic fields. A 
central closure problem in MFD theories is to compute the mean emf £ and express it in terms 
of the mean field itself. In the two-scale approach one assumes that the mean field is spatially 
smooth over scales bigger than the turbulence coherence scale /, and expresses the mean emf £ 



in terms of the mean magnetic field and its first derivative [21]. For isotropic, homogeneous, 
helical 'turbulence' in the approximation that the correlation time t is short (ideally ut/1 «: 1, 
where u is the typical turbulent velocity) one employs what is known as the First order smoothing 
approximation (FOSA) to write 



£ = u xb = a/fB-?7,VxB. (5.2) 
Here (%k = —(t/3)oj-u is the dynamo a-effect, proportional to the kinetic helicity and rj t = (r/3)u 2 /3 



is the turbulent magnetic diffusivity proportional to the kinetic energy of the turbulence [21]. 

In the context of disk galaxies, the mean velocity U is that of differential rotation. This leads 
to the fl-effect, that of shearing radial fields to produce toroidal fields, while the a-effect is crucial 
for regeneration of poloidal from toroidal fields. A physical picture of the a-effect is as follows: 
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The interstellar medium is assumed to become turbulent, due to for example the effect of super- 
novae randomly going off in different regions. In a rotating, stratified (in density and pressure) 
medium like a disk galaxy, such turbulence becomes helical. Helical turbulent motions of the gas 
perpendicular to the disk draws out the toroidal field into a loop which looks like a twisted Q.. Such 
a twisted loop is connected to a current which has a component parallel to the original toroidal 
field. If the motions have a non-zero net helicity, this parallel component of the current adds up 
coherently. A toroidal current then results from the toroidal field. Hence, poloidal fields can be 
generated from toroidal ones. (Of course microscopic diffusion is essential to make permanent 
changes in the field). This closes the toroidal-poloidal cycle and leads to exponential growth of 
the mean field. The turbulent diffusion turns out to be also essential for allowing changes in the 



mean field flux. The kinematic MFD equation 5.1, gives a mathematical foundation to the above 



picture. One finds exponentially growing solutions, of the MFD equations provided a dimension- 



less dynamo number has magnitude D = \aoGh 3 7]~ 2 \ > D cr i t ~ 6 [11, [2^ ], a condition which can be 
satisfied in disk galaxies. (Here h is the disk scale height and G the galactic shear, ao typical value 
of a, and we have defined D to be positive). The mean field grows typically on time-scales a few 
times the rotation time scales, of order 3 - 10 x 10 8 yr. 

This picture of the galactic dynamo faces several potential problems. Firstly, while the mean 
field dynamo operates to generate the large-scale field, the fluctuation dynamo is producing small- 
scale fields at a much faster rate. Also the correlation time of the turbulence measured by ut/1 is 
not likely to be small in turbulent flows. So the validity of FOSA is questionable. Indeed, based 
on specific imposed (kinematic) flow patterns it has been suggested that there is no simple rela- 
tion between a and helicity of the flow; see [51]. In order to clarify the existence of a-effect and 
turbulent diffusion as outlined above, and its R m dependence, even in the kinematic limit, we have 
recently measured £ directly in numerical simulations of isotropic, homogeneous, helical turbu- 
lence [p4j]. These simulations reach up to a modest Rm ~ 220. We find, somewhat surprisingly, 
that for isotropic homogeneous turbulence the high conductivity results obtained under FOSA are 
reasonably accurate up to the moderate values of R m that we have tested. A possible reason for this 
might be that the predictions of FOSA are very similar to a closure called the minimal t approxi- 
mation (MTA) [^7|, ||] where the approximation of neglecting nonlinear terms made in FOSA is not 
done, but replaced by a closure hypothesis. But MTA is also not well justified, although numerical 



simulations (for R m < 300) support some aspects of this closure [|28|]. Interestingly, this agreement 
of a and ij t directly measured from the simulation, with that expected under FOSA, is obtained 
even in the presence of a small-scale dynamo, where b is growing exponentially. This suggests that 
the exponentially growing part of the small-scale field does not make a contribution to the mean 
emf £, correlated with the mean field B. It is essential to extend these results to even higher values 
of R m , but these preliminary results are quite encouraging. 

Another potential problem with the mean field dynamo paradigm is that magnetic helicity 
conservation puts severe restrictions on the strength of the ar-effect The magnetic helicity 
associated with a field B = V x A is defined as Hj = /A.B dV, where A is the vector potential 



[21]. Note that this definition of helicity is only gauge invariant (and hence meaningful) if the 
domain of integration is periodic, infinite or has a boundary where the normal component of the 
field vanishes. Hj measures the linkages and twists in the magnetic field. From the induction 
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equation one can easily derive the helicity conservation equation, dHj/dt = -2j] JB • (V x B) dV. 
So in ideal MHD with rj = 0, magnetic helicity is strictly conserved, but this does not guarantee 
conservation of Hj in the limit rj — > 0. For example magnetic energy, whose Ohmic dissipation 
is governed by (dEB/dt)j ou \ e - -t](4k/c 2 ) J J 2 dV, can be dissipated at finite rates even in the limit 
77 — > 0, because small enough scales develop in the field (current sheets) where the current density 
increases with decreasing 77 as <x rf 1 ^ as r\ — > 0. Nevertheless, since helicity dissipation rate has 
a milder dependence on the current density (its rate being proportional to J • B and not J ), in 
many astrophysical conditions where R m is large (rj small), the magnetic helicity Hj, is almost 
independent of time, even when the magnetic energy is dissipated at finite rates. 

The operation of any mean-field dynamo automatically leads to the growth of linkages between 
the toroidal and poloidal mean fields and hence a mean field helicity. In order to satisfy total helicity 
conservation this implies that there must be equal and oppositely signed helicity being generated 
in the fluctuating field. What leads to this helicity transfer between scales? To understand this, 
we need to split the helicity conservation equation into evolution equations of the sub-helicities 
associated with the mean field, say H T - J A ■ B dV and the fluctuating field hj = J a ■ b dV. The 
evolution equations for Hj and hj are [Q] 



dHj_ 

dt 



r C4 n dh T f C4n 

I 2S-BdV-2?i I — J-BdV; = - I 2S-BdV-2r] I — j-bdV. 



Here, we have assumed that the surface terms can be taken to vanish (we will return to this issue 
below). We see that the turbulent emf £ transfers helicity between large and small scales; it puts 
equal and opposite amounts of helicity into the mean field and the small-scale field, conserving the 
total helicity Hj = Hj + hj. Note that in the limit when the fluctuating field has reached a stationary 
state one has dhjldt — > and f& ■ B dV = -2rj J(4n/c)j ■ b dV which tends to zero as rj — > for 
any reasonable spectrum of small scale current helicity. Therefore, the component of the turbulent 
electromotive force along the mean field is catastrophically quenched in the sense that it tends to 



zero as 77 — > 0, or for large R m , a feature also borne out in periodic box simulations [19, 25 ]. 

To make the above integral constraint into a local constraint requires one to be able to define a 
gauge invariant helicity density for at least the random small-scale field. Such a definition has in- 



deed been given, using the Gauss linking formula for helicity In physical terms, the magnetic 
helicity density h of a random small scale field (in contrast to the total helicity hj), is the density 
of correlated links of the field. This notion can be made precise and a local conservation law can 



be derived (see [26] for details) for the helicity density h, 



-^ + V-F = -2£-B-277(47r/c)j-b. (5.3) 
ot 

where F gives a flux density of helicity. (For a weakly inhomogeneous system, h is approximately 
a b in the Coulomb gauge.) In the presence of helicity fluxes, we have in the stationary limit, 
& ■ B = -(1/2)V • F - 7/j • b, and therefore £ • B need not be catastrophically quenched. Large scale 
dynamos then seem to need helicity fluxes to work efficiently [ [29| , |30| , . 

This conclusion can be understood more physically as follows: As the large-scale mean field 
grows the turbulent emf £ is transferring helicity between the small and large scale fields. The large 
scale helicity is in the links of the mean poloidal and toroidal fields, while the small scale helicity 
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is in what can be described as "twist" helicity of the small-scale field. Lorentz forces associated 
with the "twisted" small-scale field would like to untwist the field. This would lead to an effective 
magnetic a-effect which opposes the kinetic a produced by the helical turbulence (see below). The 
cancellation of the total or-effect is what leads to the catastrophic quenching of the dynamo. This 
quenching can be avoided if there is some way of transferring the twists in the small scale field out 
of the region of dynamo action, that is if there are helicity fluxes out of the system. 

Blackman and Field (in [29]) first suggested that the losses of the small-scale magnetic helicity 
through the boundaries of the dynamo region can be essential for mean-field dynamo action. Such 
a helicity flux can result from the anisotropy of the turbulence combined with large-scale velocity 
shear or the non-uniformity of the a-effect p9|]. Another type of helcity flux is simply advection 



of the small scale field and its associated helicity out of the system, with F = h\J [30]. This effect 
naturally arises in spiral galaxies where some of the gas is heated by supernova explosions produc- 
ing a hot phase that leaves the galactic disc, dragging along the small-scale part of the interstellar 
magnetic field. 

In order to examine the effect of helicity fluxes in more detail, one also needs to fold in a 
model of how the dynamo co-efficients get altered due to Lorentz forces. Closure models either 
using the EDQNM closure or the MTA or quasi-linear theory, suggest that the turbulent emf gets 
re-normalized, with a = ax + aM, where ckm = (t/3) < b.Vxb > /(4np), is proportional to the small- 
scale current helicity and is the magnetic alpha effect mentioned above [^, ^7], 32]. The turbulent 
diffusion rjj is left unchanged to the lowest order, although to the next order there arises a non- 
linear hyperdiffusive correction [16]. Some authors have argued against an % contribution, and 
suggest instead that the a-effect can be expressed exclusively in terms of the velocity field, albeit 
one which is a solution of the full momentum equation including the Lorentz force [33]. To clarify 
this issue, we have studied the nonlinear cr-effect in the limit of small /? m ,Re <sc 1 using FOSA 
applied to both induction and momentum equations [|4|]. We show explicitly in this limit that one 
can express a completely in terms of the helical properties of the velocity field as in traditional 
FOSA, or, alternatively, as the sum of two terms, a so-called kinetic a-effect and an oppositely 
signed term proportional to the helical part of the small scale magnetic field, akin to the above 
closures. 

Adopting a = czk + um as above, one can now look for a combined solution to the helicity 
conservation equation (5.3), and the mean-field dynamo equation [^, 35], after relating the current 
helicity arising in % to the magnetic helicity density h. The effect of the advective flux in resolving 
the quenching of the dynamo was worked out in detail in Ref [[30|]. In the absence of an advective 
flux, the initial growth of magnetic field is catastrophically quenched and the large-scale magnetic 
field decreases at about the same rate as it grew. The initial growth occurs while the current helicity 
builds up to cancel the kinetic a-effect. However, even a modest advective flux compensates the 
catastrophic quenching of the dynamo and the mean field stays steady at about 10% of the equipar- 
tition value. (Excessive advection, however, hinders the dynamo as it removes the mean field from 
the dynamo active region). The vertical advection of magnetic helicity by galactic fountain flow 
can therefore resolve in a straightforward fashion the catastrophic quenching of nonlinear mean- 
field galactic dynamos. Further, advection of small-scale magnetic fields may help the mean-field 
dynamo action in other ways. For example, there is a concern that the back-reaction of Lorentz 
forces due to the rapidly growing small-scale field could quench turbulent transport processes, sim- 
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ply by suppressing the required Lagrangian chaos. In this case their advection out of the galaxy 
will still allow the dynamo to operate efficiently. 

We have emphasized so far the role of helicity for the dynamo generation of large-scale fields. 
Intriguingly, several recent simulations show that the generation of large-scale fields may arise even 
in non-helical turbulence in the presence of strong enough shear [|36|]. It is at present an important 
open question as the exact cause of such large-scale field generation. 

6. Magnetizing the IGM through outflows 

We have already mentioned that outflows from supernovae and AGN could seed galaxies and 
clusters with magnetic fields. We now consider whether such processes could in fact magnetize the 
general intergalactic medium and if so at what levels. 

In this context it is interesting to note that there is growing evidence for metals in the mod- 
erately overdense intergalactic medium (IGM) up to redshifts z ~ 5 or so (cf. and references 
therein). The metals detected in the IGM can only have been synthesized by stars in galaxies, 
and galactic outflows energized by supernovae (SNe) are the primary means by which they can 
be transported from galaxies into the IGM. If galactic ISM is also magnetized due to some form 
of dynamo action, then such outflows will also transport magnetic fields into the IGM. One can 
roughly estimate the IGM field resulting from such outflows, by using magnetic flux conserva- 
tion under spherically symmetric expansion; that is B see d - (piGM/piSM^^gai- For B ga \ - 3fiG, 
and Pigm/Pism - 10" 6 , one gets B see d - 0.3 nG. In a similar fashion outflows from AGN can also 
transport magnetized plasma into the IGM. 

The level of such magnetization cannot be estimated at present with any great certainty, since 
it depends not only on a host of parameters to do with the energetics of the outflows but also on the 
evolution of the magnetic field in the outflow, whether it is amplified by dynamo action or decays 
by generating decaying MHD turbulence. Nevertheless there are several attempts [ |38| ] where some 
model of the outflows is tied in with a model of magnetic field behaviour. For example Bertone 
et al (in [|38|]) find the outflows affected regions can have fields ranging from 1(T 12 - 1(T 8 G, not 
much of a volume filling at z = 3 but a significant volume filling by z = 0. Due to outflows from 
AGN, Furlanetto and Loeb (in [0]) estimate that by a redshift z ~ 3, about 5%-20% of the IGM 
volume is filled by magnetic fields and the magnetic pressure in these regions could be comparable 
to the thermal pressure of the the photoionized IGM (at T ~ 10 4 K). 

It should also be noted that such outflow scenarios are constrained by the fact that they should 
not dynamically perturb the IGM responsible for the Lyman-a forest absorption at z = 3. In this 
context scenarios involving outflows from small mass galaxies and filling the IGM at high redshifts 
are more favored, for magnetizing the high z IGM significantly without perturbing the Lyman- 
a forest. 

Of course, as mentioned earlier, one still needs to mix the magnetized and unmagnetized media 
and what this does to both the field strength and coherence scale is as yet unclear. This would be 
relevant if we wish to use the IGM field to explain galaxy and cluster magnetism. Nevertheless, the 
magnetization of a significant volume of the IGM by z = could impact significantly on other 
astrophysical phenomena. For example significant IGM magnetic fields would perturb cosmic ray 
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propagation and is relevant to the issue of whether one can do high-energy cosmic ray astronomy 
at all. 



7. Primordial magnetic fields from the early universe? 

We have so far concentrated on the hypothesis that magnetic fields observed in galaxies and 
galaxy clusters arise due to dynamo amplification of weak seed fields. An interesting alternative 
is that the observed large-scale magnetic fields are a relic from the early Universe, arising perhaps 
during inflation or some other phase transition (p0|] and references therein). It is well known that 
scalar (density or potential) perturbations and gravitational waves (or tensor perturbations) can be 
generated during inflation. Could magnetic field perturbations also be generated? 

Indeed inflation provides several ideal conditions for the generation of a primordial field with 



large coherence scales [|40|]. First the rapid expansion in the inflationary era provides the kinemat- 
ical means to produce fields correlated on very large scales by just the exponential stretching of 
wave modes. Also vacuum fluctuations of the electromagnetic (or more correctly the hypermag- 
netic) field can be excited while a mode is within the Hubble radius and these can be transformed to 
classical fluctuations as it transits outside the Hubble radius. Finally, during inflation any existing 
charged particle densities are diluted drastically by the expansion, so that the universe is not a good 
conductor; thus magnetic flux conservation then does not exclude field generation from a zero field. 
There is however one major difficulty; since the standard electromagnetic action is conformally in- 
variant, and the universe metric is conformally flat, one can transform the evolution equation for 
the magnetic field to its flat space version. The field then field always decreases with expansion as 
I /a 2 , where a(t) is the expansion factor. 

Therefore mechanisms for magnetic field generation need to invoke the breaking of conformal 
invariance of the electromagnetic action, which change the above behaviour to B ~ 1 [a e with typi- 
cally e <sc 1 for getting a strong field. Since a{t) is almost exponentially increasing during slow roll 
inflation, the predicted field amplitude is exponentially sensitive to any changes of the parameters 
of the model which affects e. Therefore models of magnetic field generation can lead to fields as 
large as B ~ 10~ 9 G (as redshifted to the present epoch) down to fields which are much smaller than 
that required for even seeding the galactic dynamo. Note that the amplitude of scalar perturbations 
generated during inflation is also dependent on the parameters of the theory and has to be fixed 
by hand. But the sensitivity to parameters seems to be stronger for magnetic fields than for scalar 
perturbations due to the above reason. 

Another possibility is magnetic field generation in various phase transitions, like the elec- 
troweak transition or the QCD transition due to causal processes. However these generically lead 
to a correlation scale of the field smaller than the Hubble radius at that epoch. Hence very tiny 
fields on galactic scales obtain, unless helicity is also generated; in which case one can have an 
inverse cascade of energy to larger scales pi]]. 

If a primordial magnetic field with a present-day strength of even B ~ 10~ 9 G and coherent on 
Mpc scales is generated, it can strongly influence a number of astrophysical processes. For exam- 
ple, such primordial magnetic fields could induce temperature and polarization anisotropics in the 



Cosmic Microwave Background (CMB) (see [42] for a review). The signals that could be searched 



for include excess temperature anisotropics (from scalar, vortical and tensor perturbations), B- 
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mode polarization, and non-Gaussian statistics [ff3j|. A field at a few nG level produces tempera- 
ture anisotropics at the 5//K level, and B-mode polarization anisotropics 10 times smaller, and is 
therefore potentially detectable via the CMB anisotropics. An even smaller field, with Bo ~ 0.1 
nG, could lead to structure formation at high redshift z > 15, and hence naturally impact on the 
re-ionization of the Universe [44]. A 0.1 nG field in the IGM could also be sheared and amplified 
due to flux freezing, during the collapse to form a galaxy and lead to the few fiG field observed in 
disk galaxies (cf. pB]]). Ofcourse, one may still need a dynamo to maintain such a field against 



decay and/or explain the observed global structure of disk galaxy fields [22]. Weaker primordial 
fields can also provide a strong seed field for the dynamo. Overall, it is interesting to continue to 
look for evidence of such a primordial field. 



8. Conclusions 

We have presented a brief overview of issues related to the origin of cosmic magnetic fields. 
Battery mechanisms produce only a small seed field which needs to be amplified by a dynamo. 
How fluctuation and mean-field dynamos work is at present under intense scrutiny. Basic ideas 
of turbulent dynamos are in place. But a detailed understanding of the saturation of fluctuation 
dynamos, and the linear/nonlinear behaviour of turbulent transport co-efficients associated with 
MFD's, is still challenging. There is also increasing interest in finding natural mechanisms for 
primordial field generation in the early universe, and their observational consequences. Our knowl- 
edge of galactic, cluster and IGM fields has come from observing the synchrotron polarization and 
Faraday rotation in the radio wavelengths. Instruments like the SKA will expand this knowledge 
enormously p6[], and will be crucial to further probe the magnetic universe. 
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